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Abstract 
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the reduced SU{?)) flavour symmetry is explicitly broken by the general cur- 
rent quark mass term with m„ 7^ rhd 7^ rhs- 
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1. Introduction 



Phenomenological parametrizations based on some simple ansatz with solid 
symmetry grounds are frequently used in low-energy QCD. One of the most 
common and important outputs of this approach is to get a clue of how 
high energy QCD may influence low energy observables. Unfortunately, in 
spite of all remarkable successes of the QCD sum rules method or chiral 
perturbation theory (2j, this picture is still far away from being completed. 

Some features of the large distance hadron dynamics can be understood 
in the framework of effective chiral Lagrangians written in terms of quark 
degrees of freedom They are efficient for the description of spontaneous 
chiral symmetry breaking, or for the study of the quark structure of light 
mesons. The parameters of such Lagrangians can be related to the charac- 
teristics of the QCD vacuum given in form of the vacuum expectation values 
of the relevant quark bilinears or gluons (if they are included). In many re- 
spects this approach corresponds to a Landau - Ginzburg like description of 
the flavour dynamics. 

There are a number of instructive models that assume the existence of un- 
derlying multi-quark interactions and their importance for physics of hadrons. 
Well-known examples include the Nambu and Jona-Lasinio (NJL) model |3], 
where the four-fermion interactions have been used to study dynamical break- 
ing of chiral symmetry^; the instanton inspired models |H|, where 2iVj-quark 
interactions {Nf is the number of quark flavours) offer a possible framework 
to discuss the f/A(l) problem 0; the potential-type quark models which are 
successfully applied to the evaluation of hadronic parameters pUj . 

In this Letter, we propose to extend the phenomenologically interesting 
three-flavour quark model which combines the chiral U{3)l x U{3)r NJL-type 
Lagrangian with the 't Hooft six-quark determinant (NJLH), by supplying it 
with flavour mixing eight-quark interactions. The original NJLH Lagrangian 
gives a good description of the pseudoscalar nonet, especially the r] and 77' 
masses and mixing JT], and in this form the model has been widely and 
successfully explored at the mean- field level [T^-fTlj. 

This approximation was refined by works of H. Reinhardt and R. Alkofer 
[T3] . who used the functional integral method to bosonize the model. This 
approach hinges decisively on the stationary phase asymptotics of the gen- 

^Later on a modified form of this interaction lias been used to derive the QCD effective 
action at long distances El d • 
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erating functional and allows to calculate the contribution of the classical 
path already at lowest order. This lowest order result sums all tree diagrams 
in the perturbative series in powers of the coupling constant of the 't Hooft 
interaction [T^ . 

The functional treatment of the model reveals one essential problem: the 
model has actually several classical trajectories which belong to the interval 
of the functional integration, and therefore contribute to the integral 
If one takes them into account, the effective potential of the theory gets 
unbounded from below, i.e., the system does not have a ground state. 

We argue here that this drawback of the NJLH model can be removed. 
The eight-quark interactions added to the original Lagrangian reduce (under 
given conditions) the number of stationary phase trajectories to one and, as a 
result, the theory has a stable global minimum, attributed to a spontaneous 
symmetry breakdown. It should be remarked that the stationary phase equa- 
tions which appear in this approach are of cubic order and have, in general, 
more than one admissible solution. We obtain inequalities for coupling con- 
stants to distinguish those solutions further and show that these constraints 
can be finally understood as the stability criteria of the whole system. We 
consider the most general eight-quark spin zero interactions invariant under 
U{3)l X U{3)r chiral symmetry and assume that current quarks have realis- 
tic masses: m„ 7^ rh^ 7^ rfi^. It is shown that our result is independent both 
of the specific form of eight-quark interactions and values of current quark 
masses. 

Let us note that only one type of the eight-quark interactions considered 
are flavour mixing (the first part of the eight-quark Lagrangian studied here, 
i.e., £1 in eq.(jD)) and have been used previously [17j in a different context, 
namely to introduce OZI-violating effects [18j in a NJL-type model with the 
U\{1) anomaly term inspired by the works of Di Vecchia and Veneziano 
and independently by Rosenzweig, Schechter and Trahern [20] • Recently, by 
describing the properties of nuclear matter with two-flavour NJL models, 
eight-fermion interactions of the £i-type have been also analyzed in fIT\ . 

2. The model 

The dynamics of the model considered is determind by the Lagrangian den- 
sity 

^es = q{il^df, - m)q + /:njl + -Cr + (1) 
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where it is assumed that quark fields have colour {Nc = 3) and flavour {Nj = 
3) indices. The current quark mass, m, is a diagonal matrix with elements 
diag^-fhu, rhd, rhg) , which explicitly breaks the global chiral SUl{S) x SUr^S) 
symmetry of the Lagrangian. The flavour symmetry of the model becomes 
SU{3), if rhu = rha = rhs] and one gets the reduced symmetries of isospin 
and hypercharge conservation, if = rhd 7^ rhg. Putting 7^ 7^ m^, 
one obtains the most general pattern of the explicit symmetry breakdown in 
the model. 

We suppose that quark verticies are effectively local, this being a fre- 
quently used approximation. Even in this essentially simplified form the 
Lagrangian has all basic ingredients to describe the dynamical symmetry 
breaking of the hadronic vacuum and find its stability condition. 

The interaction Lagrangian of the NJLH model in the scalar and pseu- 
doscalar channels is given by two terms 



The first one is the Ul{?)) x Ur{?)) chiral symmetric interaction specifying the 
local part of the effective four-quark Lagrangian in channels with quantum 
numbers = O"*", 0~. The Gell-Mann flavour matrices A^, a = 0, 1, . . . , 8, 
are normalized such that tr(AaAb) = 26ab- The second term represents the 
't Hooft determinantal interactions [H]. The matrices Pl.r = (1 =F75)/2 are 
projectors and the determinant is over flavour indices. The determinantal 
interaction breaks explicitly the axial Ua{^) symmetry [22] and Zweig's rule. 

The new feature of the model is the inclusion of U{3)l x U{3)ji symmetric 
eight-quark forces, which we add to the standard NJLH Lagrangian to obtain 
the stable ground state. They are described by the term Csq = where 



The flavour indices i, j, . . . = 1, 2, 3 = u, d, s, and gi, g2 stand for the various 
symmetric eight-quark coupling strengths. The first term £1 coincides with 
the OZI- violating eight-quark interactions considered in ^7]. The second 
term C2 represents interactions without violation of Zweig's rule. Csq is the 
most general Lagrangian which describes the spin zero eight-quark interac- 
tions without derivatives. It is the lowest order term in number of quark 



£h = fi;(det qPiq + det qPRq). 



(2) 
(3) 



t^i = 8gi [{qiPRqm){qmPLqi)f , 

C2 = ^Qg2{qiPRqm){qmPLqj){qjPRqk){qkPLqi)- 



(4) 
(5) 
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fields which is relevant to the case. We restrict our consideration to these 
forces, because in the long wavelength limit the higher dimensional operators 
are suppressed. 

Large Nc arguments can be also used to justify this step if the dimension- 
full coupling constants [G] = M'^ [k] = M-^ [gi] = [5(2] = M"^ count at 
large A^^^ as G ~ l/N^, n ~ , gi,g2 ~ V^c • ^^is case the NJL inter- 

actions (j21) dominate over £h and Cgq at large N^, as it should be, because 
Zweig's rule is exact at Nc = 00. On the other hand, with these counting 
rules the Lagrangians £h and Csq contribute at the same Nc order, thus the 
effects coming from them are comparable and must be considered together^. 

It is clear that our considerations are also relevant if the multi-quark inter- 
actions create a hierarchy j2Sl similar to the hierarchy found within the gluon 
field correlators [21]. In this case the lowest four-quark interaction forms a 
stable vacuum corresponding to spontaneously broken chiral symmetry. The 
higher multi-quark interactions in the hierarchy must not destroy this state, 
otherwise they would be as important as the lowest order terms. Since how- 
ever the 't Hooft interaction, which is the next term in the hierarchy, destroys 
the ground state [in], one cannot truncate the tower of multi-quark interac- 
tions at this level. The next natural candidate is the eight-quark term Cgq- 
We show that its inclusion is sufficient to stabilize the ground state. 



3. The eight-quark term at work 

The many-fermion vertices of Lagrangian CeS can be presented in the bilinear 
form by introducing the functional unity J3] in the vacuum-to-vacuum am- 
plitude of the theory. The specific details of this bosonization procedure are 
given in our recent work fTHj. The new interaction term C^q, which we add 
now to the effective quark Lagrangian, does not create additional problems, 
and the method can be simply extended to the present case. This is why we 
take as a starting point the corresponding functional integral already in its 
bosonized form 

Z = VqVq'[[Vaa'[[V(f)a exp (i / d^x£g(g, g, a, 0) j 

^Let US note that our counting for g\ differs from the prescription of paper |17j where 
51 ~ 
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+ 00 

X f U^Sall Vpa exp U I d^x£,(a, 0, A; s, p)) , (6) 

where 

£^ = qii^'^d^ -m-a- i'^^(t))q, (7) 

= S„(aa + A„) + Pa(l>a + ^ (^a + Pa) 

+ y [dabedcde (SaSbScSd + 2SaSbPcPd + PaPbPcPd) 

~l~ 4 /'ace /bde ^ a SbPcPd] ■ (8) 

It is worth to observe that we did not use any approximations to obtain this 
result. 

Let us explain our notations. The bosonic fields o"a and (pa are the com- 
posite scalar and pseudoscalar nonets which will be identified later with the 
corresponding physical states. The auxiliary fields Sa and pa must be in- 
tegrated out from the effective mesonic Lagrangian C^- We assume that 
a = <7a\a, and so on for all bosonic fields a, 0, s,p. The quarks obtain their 
constituent masses m = rriaXa = diag(m^j, rrid, m^) due to dynamical chiral 
symmetry breaking in the physical vacuum state, = ma — rha. The totally 
symmetric constants Aabc are related to the flavour determinant, and equal 
to 

1 
3! 

The eight-quark interactions change drastically the semi-classical asymp- 
totics of the functional integral over Sa,Pa in &, as compared to the case, 
when gi,g2 = 0. To see this one should first find all real stationary phase 
trajectories s^* = Sa{cr, 4>), Pa = Pa{cr, 4>) given by the equations 

^ = 0, ^ = 0. (10) 

OSa dpa 

We seek these solutions in form of expansions in the external mesonic fields, 

= ha + h'abCTb + h'abcCTbac + h^abcMc + • • • 
Pa =h^aUl> + h'^l'Pb^, + --- (11) 



■^abc r)| ^Jjfcf mni ('^a)im('^6)jrn (-^c) fci • (9) 
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The coefficients h^^^ depend on the couphng constants G, k, gi, g2 and quark 
masses A^. The higher index coefficients h^^^ are recurrently expressed in 
terms of the lower ones. The one-index coefficients ha are the solutions of 
the following system of cubic equations 

A„ + Gha + — Aabchhc + — hahl + — dabedcdehhchd = 0. (12) 

The trivial solution ha = 0, corresponds to the perturbative vacuum A^ = 0. 
There are also non-trivial ones. In accordance with the pattern of explicit 
symmetry breaking the mean field A^ can have only three non-zero compo- 
nents at most with indices a — 0,3, 8. It means that in general we have a 
system of only three equations to determine haXa — diag(/iu, ha, hs) 

Ghu + A„ + ^ haK + J huihl + hl + hi) + ^hl = 0, 
Gha + Ad + ^ h^K + I ha{hl + hi + h^) + ^-^ h\ ^ 0, (13) 
Ghs + A. + ^ h^ha +^-^K{hl + hl + hl) + ^-^hl = Q. 

Our aim now is to show that parameters can be fixed in such a way that 
this system will have only one real solution. We start by summing the first 
two equations, which leads to the cubic equation 

+ tx^ b, 

t = (SG +'^hs + y\gi + 8^2) + ^g^hl) , 

91+92 \ 2 J 

^^_8(A^^ (14) 
91 + 92 ^ ^ 

where x = hu + ha, y = hu - ha- 

Note that deviations of the variable y from zero are a measure of isospin 
breaking effects due to electromagentic forces, as the difference hu — hd does 
not vanish for m„ 7^ rhd- The function t{y, hg) has a minimum (if gi > and 
gi + 3(72 > 0) at y = and hg = —k/ (Sgi), thus the inequality t > always 
holds for coupling constants fixed by 
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In this case the cubic equation has for any given value of b just one real root 

Since 6 < (provided that A„ + A^^ > 0), this function is negative. Its 
minimum is located at the point y = 0, hg = —K/{8gi), and the surface 
X = is an asymptotic one to 

Subtracting the second equation from the first one we obtain a quadratic 
equation with respect to x. Its solutions are given by 



X(2) = ±. 



(17) 

Since we only allow real solutions, the following inequality must hold 



%G-'^hs + 2g,hl + y\g, +g2)<- (A, - AJ . (18) 

y 

For definiteness, we suppose that A^ — A„ > 0. This assumption represents 
one of two possible alternatives. Our final mathematical conclusions do not 
depend on the choice made. However, it is not obvious which of them should 
be required physically. Next, the function fihg) = 8G — Khs/2 + 2gih1 > 0, 
since the minimum value fmin = ~ ^ I (32(?i) > in the parameter region 
of (fT^ . Therefore y ranges over the half-open interval < y < ymax{hs). 
The lower bound y = is an asymptotic surface for the function X(2). The 
upper bound ymax{hs) is the unique real solution of the equation y^{gi + 
92) + yf{hs) + 8(A^, - Ad) = 0. It follows that ymax oc (A^ - A„), i.e., 
the electromagnetic forces which are responsible for the isospin symmetry 
breaking determine the length of the segment [0,?/max], which is relatively 
small as compared with intervals determined by the strong interaction. As a 
consequence a negative branch of the function X(2) grows rapidly with y from 
—00 at y = up to at ?/ = Vmax- On the contrary, functions a;(i)(y) and 
X(3)(?/) (see eq. (fT^ below) remain almost unchanged in the interval < y < 
ymaxi because here the strong driving forces totally cover electromagnetic 
effects. 

Let us consider now the third equation which yields 



vjhs^y) 
^ -{k + 8gihs 



-(3) = ±4 , ^^^^^^^ , (19) 
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Figure 1: The label x stands for the curves plotted as function of y at 

fixed hs = -0.03 GeV^ The bellshaped curve corresponds to the two branches of 
X{2), its peak is located at ymax{hs) = 8.8 • 10~''GeV^. The thick line indicates 
the negative branch of curve X(3) and also the curve X(^i^, which are degenerate at 
this scale, the upper line is the positive branch of X(3). The two branches of X(3) 
meet at y ~ 0.14 GeV^, outside the indicated range in the plot. The real solution 
of the system (|13|) is indicated by the point P. The corresponding parameters are 
G = 6.1GeV^2^ 91=92 = 4•103GeV-^ k = -705GeV-^ A„ = 313 MeV, = 
318 MeV, A, = 345 MeV. 

where we have introduced the notation 

v{hs, y) = {gi + 2g2)hl + ^{8G + g^y'') - ^y" + 4A,. (20) 

2 lb 

The expression under the square root is positive, if conditions 

v{K,y)>^, K + 8^i/i, <0, (21) 

are fulfilled. The alternative case does not have solutions, since we assume 
that As > and k < (phenomenological requirements). Inequalities (|2H) 
hold with hs belonging to the half-open interval /i™™ ^ hg < h^°'^ . Here 
j^max _ _f^^(^gg_^^ > g^j^j j^rnm ^ g_ rpj^g lower bouud is a solution of the 

equation v{hs, y) = 0. This cubic equation has only one real root which is 
negative for 

gi + 2(72 > 0, 8G + giy^ > 0, 4A, - ^ y^ > 0. (22) 



9 



0.4 



0.2 

X 



-0.2- \ I 

-0.4i I .1.1 - 

-0.1 -0.05 0.05 0.1 

hs 

Figure 2: The curves X(i) of eq. (fT6|) (small dashed line) and X(3) of eq. (|19j) (solid 
line) are shown as functions of hg for the parameter set of fig.l, at fixed y = 
4.2 • lO-^GeV^. The solution P has /i„ = -0.02747, ha = -0.02789, = -0.03 
in units GeV^. The vertical dashed line corresponds to hg = —K/{8gi). 

Under the assumptions made above these inequalities are obviously fulfilled. 

We illustrate the case in two figures. The y-dependence is shown in fig.l. 
Since X(2) is a monotonic function of y in the region x < 0, < y < ymax at 
any fixed value of hg, the question whether the system (jl3|) has one or more 
solutions is now reduced to a careful check of the number of intersections for 
curves and X(3) as functions of hg at a fixed value of y. Actually, for 
this purpose one can choose any value of y from the interval < ?/ < ymax-i 
because functions and X(3) are almost insensitive to this value. 

In fig. 2 we show X(^i){hs,y) and X(^3){hs,y) as functions of hg, at fixed y 
given by the solution P of fig.l. It is quite easy to verify that the line hg = 
f^max^ being the asymptote for the curve given by eq. lfT^ . crosses the other 
curve p6|) in its minimum, dividing it in two monotonic parts. Thus, both 
functions decrease monotonically with increasing hg in the third quadrant of 
the Cartesian coordinates system formed by the line hg = h'^"'^ and the axis 
of abscissas. The curves have only one intersection, which corresponds to an 
unique solution of the system (fT^. 

The fact that the cubic equations ()13p have only one set of real roots 
over a certain range of values of parameters G, k, (71, (72, is crucial for the 
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ground state of the theory: it makes the vacuum globally stable^. 

Unfortunately, we merely can find the solution /i^, /i^) numerically, 
apart from the simplest case with the octet flavour symmetry, where current 
quarks have equal masses m„ = rhd = rhg, and the system (fT^ reduces to a 
cubic equation for only one variable hu = = hg 



K + 



K 2 4G 4A 
m ""^ 3A 3A 







with X = Qi + (2/3)(72- Making the replacement 
obtains from 



hi + t'K = b', 



(23) 

hu — k/(36A), one 
(24) 



where 



t' 



G 



-)' 

24 A 7 



b' 



A { K 

3 jsGA 



G 



2 / K 



A' 



3 V24A. 

It is clear now that this cubic equation has one real root, if t' > 0, i.e.. 



G 

A> 



K 



24A 



(25) 



(26) 



In this particular case the proof of existence and uniqueness of the solution 
is straightforward. Let us also note that the solution found above for the 
general case deviates not much from the case with octet symmetry, i.e., we 
have approximately hu — hd — hg. 



4. Effective potential 

Since the system of equations (fTUj) can be solved, we are able to obtain 
the semi-classical asymptotics of the integral over Sa,Pa in One has 

the following result which is valid at lowest order of the stationary phase 
approximation 

Z[a,(f),A] = / JIX'SaH^Pa exp U / d'^xCr{a,(f),A;s,p)] 

^ a. o. ^ •' ^ 

^Let us recall that putting gi = g2 = 0, one obtains from (|13|) the system of quadratic 
equations to find hu,hd,hs. It has been shown in 16 that such equations have two real 
solutions (for a physical set of parameters) in the SU (3) case and three real solutions in 
the SU{2) X U{1) case. This is exactly the underlying reason for the vacuum instability. 
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j=l V / 

where n is the number of real solutions {s^a,p^^)^^^ of eq. (jl(J|) . 

The information about the vacuum state is contained in the effective 
potential of the theory. To obtain it let us consider the linear term in the cxa 
field. The resulting contribution, as it follows from eq. (|?fj) . is 

Z ~ exp (^z J d'x X: /i|fV, + . . . j (28) 

This part of the Lagrangian is responsible for the dynamical symmetry break- 
ing in the multi-quark system and taken together with the corresponding part 
from the Gaussian integration over quark fields in eq.lP) leads us to the gap 
equations (for each of quark's flavours i = u,d, s), 

$:# + -^m,Jo(mn=0, (29) 
i=i 

where Jo(m^) is the tadpole quark loop contribution with a high- momentum 
cutoff A ^ 

Jo(m2) = A'' -m^,\n{l + ^] . (30) 
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Using standard techniques [23 , we obtain from the gap-equations the 
effective potential U{mi) as a function of the constituent quark masses rrii 
which corresponds, in general, to the case with n real roots. Here it is more 
convenient to use (/i„, hd, hg) as independent variables, with masses being 
determined by eqs. ()13j) . In particular, if the parameters of the model are 
fixed in such a way that eqs. p3|l have only one real solution, the effective 
potential (up to an unessential constant, which is omitted here) is 

U{K, hd, K) = ^ [AGhl + KKhdK + ^ {hlf + 3(72/^^) 

- \ (31) 
where = hi + h^, + hi hf = hf, + h^, + ht and 



vim-) = 



7r2 



(32) 
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In the specific and limited case where one deals with the octet SU (3) sym- 
metric model and = the effective potential U (m) is an even function of 
m for /t = and its plot has the standard form of the double well ( "mexican 
hat") with two symmetric minima, at m = ±mmin, and one local maximum, 
at m = 0. The 't Hooft interaction («; 7^ 0) makes this curve asymmetric: if 
K < 0, the minimum located at positive values of m gets deeper as compared 
with the other minimum at negative m, becoming therefore the global min- 
imum for the whole effective potential. It corresponds to the stable ground 
state of the system with spontaneously broken chiral symmetry. 

To appreciate the correlation found between the number of critical points 
and stability let us consider the same SU{3) symmetric model in the range 
with three real roots. In this case 



3 



Y.HI> = ~ (33) 



and we find 



mVo(m2)+AMn(l + — 



(34) 



As opposed to U (m) the potential with three real roots, described by the 
function V{m) has at most a metastable vacuum, for k/X > 0. li k/\ < 0, 
the effective potential does not have extrema in the region m > 0. In both 
cases the theory related with ¥{171} is unbounded from below and is physically 
nonsensical. 

Several special properties of the eight-quark interactions are expressed in 
these results. 

Firstly, the couplings of the considered model can always be chosen to 
fulfill the inequalities required, i.e., 

^i>0, ^i + 3^2>0, Ggi >{K/16f (35) 

and therefore the model can be simply driven to its stable regime where only 
one critical point determines the asymptotical dynamics of the system. The 
NJLH model does not have a proper mechanism for that. 

Secondly, the eight-quark terms dominate at large values of hi, making 
U{hi) positive in all directions hi ±00. As a result, the function U{hi) is 
bounded from below, and exhibits a global ground state. 
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Thirdly, the effective potential coincides at gi, g2 = with the poten- 
tial obtained in the framework of the NJLH model by the mean-field method 
[T^ . This is probably at the heart of the success of NJHL: although the limit 
91, 92 ^ in U {hi) is formally not allowed as soon as inequality (fT^ does not 
hold (instead the system is then described by the unstable potential V{hi)), 
the eight-quark terms in eq. (j3ip are not likely to destroy the results of the 
mean field approach. Nevertheless one should expect some new noticeable 
effects from it. 

We must conclude that the eight-quark interactions play a fundamental 
role in the formation of the stable ground state for the unstable system 
described by the NJLH Lagrangian. We consider this finding as the main 
result of our study. 

5. Summary and discussion 

Let us summarize what we have found. 

(1) An eight-quark extension of the conventional three-flavour NJL model 
with the explicit Ua{1) breaking by the 't Hooft determinant has been sug- 
gested. We have taken the eight quark interactions in its most general form 
for spin zero states. Eight-quark interactions prove to be essential in sta- 
bilizing the vacuum of the theory: the quark model considered follows the 
general trend of spontaneous breakdown of chiral symmetry, and possesses 
a globally stable ground state, when relevant inequalities in terms of the 
coupling constants hold. 

(2) The so ensured stability of the ground state is crucial for applications 
of the model to the study of cases in which corrections (radiative, temper- 
ature, density, and so on effects) may qualitatively change the structure of 
the theory, e.g., by turning minima in the effective potential into maxima. 
Presently the U{3)l x f/(3)/j chiral symmetric NJL model with the six-quark 
't Hooft interactions is frequently used for that. The eight-quark extension 
of the model considered here is needed for well-founded calculations in this 
field. 

(3) The eight-quark interactions are an additional (to the 't Hooft de- 
terminant) source of OZI-violating effects. They are of the same order, for 
gi ~ It is important to take them into account from the phenomeno- 
logical point of view: the details of OZI-violation are still a puzzle of non- 
perturbative QCD ^Bj . 
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